Mon. Not. R. Astron. Soc. 000, [l]-?? (2002) Printed 20 December 2012 (MN M?eX style file v2. 2) 



Electron Temperatures and Free-Electron Energy 
Distributions of Nebulae from C II Dielectronic 
Recombination Lines 

Peter J. Storey^* and Taha Sochi^^f 

^ University College London, Department of Physics and Astronomy, Gower Street, London, WCIE 6BT 



Accepted XXX. Received XXX; in original form XXX 



ABSTRACT 

A recently generated theoretical line list of C ii dielectronic recombination lines to- 
gether with observational data gathered from the literature is used to investigate the 
electron temperature in a range of astronomical objects, mainly planetary nebulae. 
The electron temperature is obtained by a least-squares optimisation using all the re- 
liable observed lines in each object. In addition, the subset of lines arising directly from 
autoionising states is used to directly determine the free-electron energy distribution 
which is then compared with various theoretical possibilities. The method described 
here can potentially determine whether there are departures from Maxwell-Boltzmann 
distributions in some nebulae, as has been recently proposed. Using published obser- 
vations of the three planetary nebulae where the relevant lines are recorded, we find 
that the data are best matched by Maxwell-Boltzmann distributions but that the 
uncertainties are sufficiently large at present that K-distributions or two-component 
nebular models are not excluded. 

Key words: planetary nebulae: general - atomic processes - methods: numerical - 
radiation mechanisms: general - ISM: abundances - stars: kinematics and dynamics. 



1 INTRODUCTION 

Recombination plays an essential role in the physical pro- 
cesses that occur in nebulae, the principal electron-ion re- 
combination processes being Radiative Recombination (RR) 
and Dielectronic Recombination (DR). Recombination of an 
electron and ion may take place through a background con- 
tinuum known as radiative recombination, or through a res- 
onant recombination process involving doubly-excited states 
known as dielectronic recombination. The latter can lead ei- 
ther to autoionisation, which is a radiationless transition to 
a lower state with the ejection of a free electron, or to stabil- 
isation by radiative decay to a lower bound state, possibly 
the ground state, with a bound electron. The RR and DR 
processes are closely linked and the difference between them 
may therefore be described as artificial; quantum mechani- 
cally they are indistinguishable. 

In section[2]of this paper we determine an electron tem- 
perature from dielectronic recombination lines of C*" in the 
spectra of a number of astronomical objects, mainly plane- 
tary nebulae, using a least squares optimisation method with 
theoretical data obtained from the recently-computed theo- 



retical line list, SSI, of Sochi & Storey (2012 \ and astronom- 



ical data gathered from the literature. The theoretical list 
was generated u sing the R-ma trix ( Berrington et tt]||1995[). 



Autostructur^ ( [Eissner et al 1974 Nussbaumer & Storey 
1978) and Emissivity (Sochi 20101 codes with an interme- 
diate coupling scheme where the lines are produced by DR 
processes originating from low-lying autoionising states with 
subsequent cascade decays. The method of formation of DR 
lines is particularly simple, often only requiring the radia- 
tive probability of the transition in question and experimen- 
tally known properties, the energy and statistical weight of 
the autoionising upper state. Even in more complex cases, 
only autoionisation probabilities are required in addition. 
This can be compared with the complex recombination and 
collisional-radiative processes involved in obtaining effective 
recombination coefficients for transitions between low-lying 



ionic states (e.g. Davey et al (2000l). Determining electron 
temperatures from DR lines can therefore provide valuable 
evidence about the temperature structure of photoionised 
nebulae. 

There is a long-standing puzzle in the physics of plan- 
etary nebulae involving the discrepancy between electron 
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temperatures and ionic abundances derived from optical re- 
combination lines (ORLs) and coUisionally excited forbid- 
den lines (CELs). Although the forbidden lines are much 
stronger than the recombination lines, they are highly de- 
pendent on temperature, making abundance determinations 
potentially uncertain. On the other hand, the recombination 
lines are weak and prone to blending and can be easily con- 
taminated by radiation from other excitation processes such 
as fluorescence. Despite all these differences, there is a com- 
mon feature between the results obtained from these lines; 
that is for all the atomic species investigated so far (mainly 
from the second row of the periodic table such as C, N, O 
and Ne) the forbidden lines in planetary nebulae normally 
produce lower ionic abundances than the corresponding val- 
ues obtained from the recombination lines. The ratio of the 
ORL to the CEL abundances, the so-called abundance dis- 
crepancy factor or ADF, is case dependent and can be a 
factor of 30 or even more. This has cast some doubt on 
the validity of the previously accepted CELs analysis re- 
sults, although the stability of CEL abundances and the 
wide variations in ORL abundances between objects sug- 
gest that we should seek the solution to the problem in the 
physics or origin of the ORLs. The abundance problem ap- 
pears to be correlated to the differences between the temper- 
atures obtained from the Balmer jump of H I and that from 
the coUisionally-excited forbidden lines where the latter is 
systematically higher than the former ( |Kholtygin| 1998| |Liu| 
2002 Tsamis et al||20d7 l. In fact, obtaining higher electron 
temperatures from forbidden lines than those deduced from 
recombination lines is a general trend in nebular studies. 

Several explanations have been proposed to justify these 
discrepancies individually or collectively, though no one 
seems to be satisfactory or universally accepted. One ex- 
planation is the sensitivity of the collisionally-excited lines 
to temperature and structure fluctuations where these fluc- 
tuations within the nebular structure result in systematic 
underestimation of the heavy element abundances deduced 
from the forbidden lines. The existence of knots depleted of 
hydrogen with high heavy element contents within the neb- 
ular gas has been proposed as a possible reason for these 
fluctuations and subsequent implications. The temperature 
inside these knots of high metallicity, and hence high opac- 
ity to stellar ultraviolet emissions and large cooling rates, is 
expected to be too low for efficient production of forbidden 
lines though it is still sufficiently high for the emission of 
recombination lines. Consequently, the recombination and 
collisional lines originate in different regions of the nebular 
gas with different elemental abundances and different tem- 
peratures. However, the presence of such knots in most or all 
planetary nebulae, as would be required to explain the sys- 
tematic nature of the observations, has not been confirmed 
observationally ( Liu et aZ|1995 Garnett fc Dinerstein|20M" 



Tsamis et aZ|2003[|Liu et ap004\ 



In a recent paper by NichoUs et al { 2012 \ it is suggested 
that this long-standing problem in planetary nebulae and 
H II regions could arise from the departure of the electron 
energy distribution from the presumed Maxwell-Boltzmann 
equilibrium condition, and that it can be resolved by as- 
suming a K-distribution for the electron energy following a 
method used in solar data analysis. The electron energy dis- 
tribution will be the subject of section [4] where we consider 
only lines originating from resonance states, that is free-free 



(FF) and free-bound (FB) transitions. We obtain a direct 
sampling of the electron energy distribution from the ob- 
servational de-reddened flux and the theoretically-obtained 
parameters such as the departure coefflcients of the involved 
autoionising states and the radiative probabilities of these 
transitions. 



2 METHOD OF ELECTRON TEMPERATURE 
INVESTIGATION 

2.1 Theory 

As indicated already, the theoretical data of the C ll di- 
electronic recombination transitions and subsequent cascade 
decay are obtained from the SSI line list of [Sochi fc Storey] 
( |2012[ ) which consists of 6187 optically-allowed transitions 
with their associated data such as emissivity and effective 
recombination coefficients. The autoionising states involved 
in the transitions of this list consist of 64 resonances be- 
longing to 11 symmetries (J — 1, 3, 5, 7, 9, 11 half even and 
J = 1, 3, 5, 7, 9 half odd) which are all the resonances above 
the threshold of C^"*" ls'^2s^ ^S° with a principal quantum 
number n < 5 for the combined electron. These include 
61 theoretically-found resonances by the K-matrix method 
( [Sochi fc Storey|2012p plus 3 experimental ones which could 
not be found due to their very narrow width. The bound 
states involved in these transitions comprise 150 energy lev- 
els belonging to 11 symmetries (J — 1,3,5,7,9 half even 
and J = 1,3,5,7,9,11 hah odd). These include 142 theo- 
retically found by R-matrix, which are all the bound states 
with effective quantum number between 0.1-13 for the outer 
electron and ^ Z ^ 5, plus 8 experimental top states which 
are the levels of the ls^2s2p(^P°)3d ^F" and *D° terms. 

The theoretical and computational backgrounds for the 
atomic transition calculations including the emissivity ther- 



modynamic model are given in Sochi (20101 and Sochi fc 
[Storey[ ( [2012[ ). The calculations were performed using an 
elaborate C ionic target in the intermediate coupling 
scheme. The list has also been validated by various tests 
including comparison to literature data related to autoion- 
isation and radiative transition probabilities and effective 
dielectronic recombination coefficients. The theoretical pa- 
rameters for the bound and resonance states were also com- 
pared to the available experimental data from the National 
Institute of Standards and Technologjj^ and found to agree 
very well both in energy levels and in fine structure splitting. 

Processes other than dielectronic recombination have 
not been considered in the atomic scattering and transition 
model of the SSI list, so the results of SSI are incomplete 
for states likely to be populated by radiative recombination 
or collisional excitation and de-excitation. Any transition 
in which the upper state has an excited ion core (usually 
2s2p(^P°)) will have negligible population by radiative re- 
combination in typical nebular conditions and is therefore 
well represented by only dielectronic recombination and sub- 
sequent cascade processes. This includes all free-free and 
free-bound transitions plus those bound-bound transitions 
involving excited ion core states. 

However, for the determination of temperature we also 
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include the 4f-3d transition, A4267 A, which is the strongest 
optical C 11 recombination line and is populated almost ex- 
clusively by radiative recombination. As discussed above, 
the formation mechanism of the dielectronic lines is very 
simple while A4267 lies at the bottom of a complex cascade 
process and our analysis provides a way, in some cases at 
least, of testing and validating the results obtained from 
A4267 A and other low-lying transitions. For example, the 
possibility of the existence of some unknown mechanism that 
overpopulates the levels of the upper state of A4267 transi- 
tion causing the enhancement of the ORL abundance may 
be ruled out if the results with A4267 are consistent with 
those obtained without A4267. 

For A4267, recombination coefBcients were taken from 



case B of Davey et al (20001. These were obtained within 



a more comprehensive theory that includes radiative and 
dielectronic recombination and all relevant collisional pro- 
cesses. They were obtained in LS-coupling rather than in- 
termediate coupling but the coupling scheme should not sig- 
nificantly affect this transition. 

2.2 Observational data 

We carried out a search for C ll recombination line data in 
the literature in which over 140 data sets related mainly to 
planetary nebulae were catalogued. All data sets that com- 
prise only bound-bound transitions with no doubly-excited 
core upper state were removed. The remaining data sets 
were subjected to a refinement process in which the flux of 
all the observational lines in each data set were normalised 
to the flux of a reference observational line in the set, which 
is usually chosen as the brightest and most reliable, while 
the emissivity of all the theoretical lines in the set were nor- 
malised to the emissivity of the corresponding theoretical 
line. The ratio of the normalised observed flux to the nor- 
malised theoretical emissivity of each line were then plotted 
on common graphs as a function of temperature on log-linear 
scales. A sample of these graphs is presented in Figure |4]for 
the planetary nebula NGC 5315. All lines that did not ap- 
proach the ratio of the reference line within an arbitrarily- 
chosen factor of 3 were eliminated. The arbitrary factor of 
3 was chosen as an appropriate limit considering practical 
factors that contribute to errors in the collection of observa- 
tional data. The refinement process also involved the utilisa- 
tion of graphs in which the ratio of theoretical emissivity to 
observational flux of all lines in a certain data set was plot- 
ted on a single graph as a function of electron temperature. 
A sample of these graphs is shown in Figure [l] 

Some lines were also eliminated for various reasons re- 
lated mainly to an established or suspected misidentification 
of the line or its intensity. For example, the wavelength of the 
alleged C ll line may not match with any known theoretical 
transition. Also, the absence of a strong line in the observa- 
tional data associated with the presence of a much weaker 
line with no obvious reason casts doubt on the identification 
of the weaker line. The line may also be eliminated because 
its intensity ratio relative to another well-established line 
does not comply with the ratio obtained from theory. A 
very few lines were also out of the wavelength range of our 
line list and hence were eliminated due to lack of theoretical 
emissivity data. 

The selected refined data sets were then subjected to a 



least squares optimisation procedure which is outlined in the 
following section. It should be remarked that the observed 
flux used in the least squares procedure is the de-reddened 
flux obtained by correcting for extinction and other sources 
of error as stated by the data source and not the raw flux 
data. Therefore, there should be no ambiguity when we use 
'observed' flux in the following sections. 



2.3 Least Squares Minimisation 

In our least squares calculations we use a single fltting pa- 
rameter, which is the electron temperature, while the ob- 
servations are the flux data of the C ll recombination lines 
that we obtained from the literature. All blended C ll lines 
in the observational list are combined by considering them 
as a single line with a single flux, while C ll lines blended 
with non-C ll lines are eliminated. To compare the theo- 
retical emissivity to the observational flux, the theoretical 
emissivity of each line is normalised to the total theoretical 
emissivity of all the lines involved in the least squares proce- 
dure, while the observational flux of that line is normalised 
to the total observational flux of these lines. The normalised 
theoretical emissivities corresponding to a particular obser- 
vational flux are added when the observational flux is given 
for a whole multiplet. 

The defined by the following equation 



JV 

= 2 



{ir - 



=1 



(1) 



is computed, where i is an index running over all the TV 
lines involved in the least squares procedure, /"° and e"' are 
the normalised observational fiux and normalised theoretical 
emissivity of line i respectively, rj is the number of degrees 
of freedom, and ajno is the variance of 7"°. This variance 



is computed from formulae given in Sochi (20121. For the 



data sets with given observational errors the reported errors 
were used while for the data sets with no reported error 
a Poisson distribution was assumed and the error on the 
observed flux was assumed to be proportional to the square 
root of the flux. In some data sets, the observational error 
was given for some lines only, and hence the average of the 
given errors was assigned to the missing ones. In some cases 
when the reported error was unrealistically small resulting 
in large , the curve was scaled to unity at the minimum 
to obtain a more realistic error estimate. 

The temperature of the object is then identifled from 
its value at the minimum , while the confldence interval is 
identifled from the values of the temperature corresponding 
to the values of Xmin ± 1 on the lower and upper sides. In 
some cases, the x'^ curve was too shallow on one side and 
hence it resulted in a broad confldence interval on that side. 

In the following section we present results for those ob- 
jects where there were sufficient adequate observations to 
derive a temperature. 



3 DERIVED TEMPERATURES 

In this section we present the astronomical objects that have 
been investigated. The objects are mainly planetary nebu- 
lae and the physical parameter of interest is the electron 
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Figure 1. The ratio of theoretical cmissivity to observational flux 
as a function of temperature for the selected C ll lines from the 
NGC 7009 spectra of [Fang Liu| ( poTTl l . 



temperature of the line emitting regions. We also include 
three objects which are not planetary nebulae for compar- 
ison, where similar techniques have been used in the past. 
The theoretical and observational data for the transitions 
used in this investigation are given in Tables [1] and |2] 



3.1 NGC 7009 

NGC 7009 is a bright planetary nebula which has a double- 
ringed complex spatial structure with a rich recombination 
spectrum and a relatively large ADF. The observational data 



for this object were obtained from Fang & Liu (2011 1 where 



9 lines, listed in Table [2] were chosen following the selection 
process. In Figure [l] the ratio of theoretical emissivity to 
observational flux is plotted against electron temperature on 
a linear-linear graph for these lines. If there were no errors in 
the observational fluxes or theoretical emissivities and the 
nebula was at a single uniform temperature, these curves 
would all intersect at the same temperature. 

The graphs, with and without A4267, are given in 
Figures |2]and[3] As seen, The first indicates a tempera- 
ture of about 5800 K while the second a temperature of 
about 5500 K, which are in good agreement. We note also 
that in Figure [l]the curves for A4267 and the DR doublet 
AA8794, 8800 A show very good agreement for T > 8000K. 
This indicates both that the theoretical emissivities of these 
lines are entirely consistent and also that the observational 
data for this object for these three lines are accurate. 

In Table [3] we list electron temperatures from the cited 
literature derived from a range of ions and three types of 
spectral features: CELs, ORLs and the Balmer and Paschen 
discontinuities. As can be seen, the values that we obtained 
here are in a broad agreement with the temperature ob- 
tained from several recombination lines. The lower electron 
temperature from optical recombination lines compared to 
the values obtained from coUisionally-excited lines is consis- 
tent with the trend of the discrepancy between the abun- 
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Figure 2. Temperature dependence of ior NGC 7009 with 
the inclusion of hne A4267, where T = 5800 K at x^in with a 
confidence interval between 4961 - 6318 K. 
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Figure 3. Temperature dependence of x^ for NGC 7009 with 
the exclusion of line A4267, where T = 5500 K at X^in ™ith a 
confidence interval between 4697 - 6156 K. 



dance and temperature results of ORL and CEL (higher 
ORL abundance and lower electron temperature as com- 
pared to the CEL abundance and temperature). If this dis- 
crepancy is caused by departures from a Maxwell-Boltzmann 
distribution by the free electrons, we might expect objects 
with the largest ADF to show the largest departures, so 
we also list ADF values from the literature in Table [S] 
NGC 7009 has the largest ADF of the objects in our sample 
of planetary nebulae. We discuss its electron energy distri- 
bution in Sec [4] 



3.2 NGC 5315 

NGC 5315 is a young dense planetary nebula in the southern 
constellation Circinus located at a distance of about 2.6 kpc 



© 2002 RAS, MNRAS 000, [T]-?? 



Electron Temperatures and Electron Distributions of Nebulae from C ll Lines 5 



4267.15 
5145.16 
6791 .47 
6800.68 



BB 
BB 
BB 
BB 



3d2D 16<-4f2F*20 

3s'4P* 6 <- 3p'4P 6 
3s'4P* 4 <- 3p'4D 4 
3s'4P* 6 <- 3p'4D 6 




Temperature(K) 



X 10 



Figure 4. Ratio of normalised observed flux to normalised theo- 
retical emissivity (NFER) versus temperature on log-linear scales 
for the planetary nebula NGC 5315. 



with an interesting complex flower shape appearance. The 
observational data, which consist of 4 selected lines indi- 
cated in Table [2] were obtained from Peimbert et al (20041. 
Figure|4]displays the ratio of the normalised observed flux to 
the normalised theoretical emissivity versus electron temper- 
ature for these lines. Our calculations, with and without 
line A4267, indicate T ~ 7400 K and T ~ 6500 K respec- 
tively. Table [3] presents some values of the electron tem- 
perature of NGC 5315 as reported in the cited literature. 
Again, the temperature is lower than that derived from the 
collisionally-excited lines. 



3.3 NGC 7027 

NGC 7027 is a compact, bright, young, high excitation plan- 
etary nebula with one of the hottest central stars known for 
a PN. Two observational data sets related to NGC 7027 



were investigated: one obtained from Baluteau et al (19951, 



and the other obtained from Zhang et al (20051. Eight lines, 



listed in Table [2) were selected from Baluteau et al. The least 
squares optimisation of these lines indicate T ~ 11100 K, in 
good agreement with values obtained by other researchers; 
a sample of which is presented in Table [3] In the case of 
Zhang et al, 20 lines, given in Table [2] were chosen. The re- 
sults of least squares, with and without line A4267, indicate 
T ~ 12500 K and r ~ 12000 K respectively. 



3.4 IC 418 

IC 418 is a bright, young, carbon-enhanced, low-excitation, 
highly symmetric planetary nebula with apparent ring struc- 
ture located at a distance of about 0.6 kpc in the constel- 
lation Lepus. Our observational data on IC 418 come from 
Sharpee et al ( 2003 1 where 22 lines, presented in Table [2j 
were selected for least squares minimisation. The plots, 
with and without line A4267, minimise at T ~ 8700 K and 
T ~ 7700 K respectively. 



3.5 NGC 2867 

NGC 2867 is a compact planetary nebula with compara- 
tively small size and fairly strong surface brightness located 
at a distance of about 0.5 kpc in the southern constellation 
Carina. The observational data of this object were obtained 
from Garcia- Rojas ei"a/| ( |2009[ ) where two knots have been 
studied: one labelled NGC 2867-1 and the other NGC 2867- 
2. Two lines, given in Table [2] were selected for the optimisa- 
tion process. The NGC 2867-1 result indicates T ~ 14300 K 
while NGC 2867-2 result indicates T ~ 16000 K. The differ- 
ence in temperature value may be caused by the diflerence 
in the physical conditions of the two knots. Table [S] presents 
electron temperatures derived in previous works from transi- 
tions of different species. As seen, our values are significantly 
higher than most of the values reported in the literature. 
However, this may be explained by the complex structure of 
this nebula and the possibility of different lines being orig- 
inating from different regions with very different physical 
conditions. 



3.6 DQ Herculis 1934 

This is a peculiar old classical galactic nova originating from 
an accreting cataclysmic variable binary system which ap- 
parently consists of a white and a red dwarf. Our obser- 
vational data of this object were obtained from |Ferland et\ 
( |1984[ ) where 2 lines, given m Table [2) were chosen for 
least squares. The result indicates T ~ 1600 K, in very good 
agreement with some values reported in the literature (refer 
to Table m notably those of |Smits| ( [l99T| ) and |Davey| ( |1995 1 
which are also derived from C 11 recombination lines. The 
higher temperatures in the literature may belong to the hot 
inner disk region, rather than the cool outer shell, where 
much higher temperatures have been derived. 

It should be remarked that the use of A 1335 A line 
is a second exception (the first is A4267 A) to our rule 
of using the BB transitions only if the upper state has a 
doubly-excited core. The upper state of the A1335 transition 
is ls'^2s2p'^ which is connected to the C'^ 



uum and ls^2s^n/ Rydberg states by two-electron radiative 
processes which are usually very weak. Hence excitation of 
A1335 A by direct radiative recombination can be neglected. 



3.7 CPD - 56°8032 

CPD - 56° 8032 is a cool late- type Wolf-Rayet star that is 
usually classified as WCIO or WCll. The star, which is lo- 
cated at about 1.3-1.5 kpc, is surrounded by a young plane- 
tary nebula with complex visible structure. Here, we try to 
infer the electron temperature of the stellar nebular wind 
surrounding CPD - 56°8032. The observational data of this 



object were obtained from De Marco et al (19971 where 13 
lines, given in Table [2] were extracted following a selection 
process. The plot indicates T ~ 17300 K. This agrees, 
within the reported error bars, with the temperature of De 
Marco et al flDe Marco||1996| IDe Marco et aZ||1996|) who 



deduced a value of 18500±1500 K for this object using a 
similar least squares approach. Our value also agrees rea- 
sonably with some of the values reported in the literature; 
a sample of which is given in Table [3] 
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3.8 He 2-113 

He 2-113 is another late- type WCIO Wolf-Rayet star sur- 
rounded by a planetary nebula with an apparent ring struc- 
ture. There are many physical similarities between CPD - 
56°8032 and He 2-113 such as age, flux and distance. These 
similarities are reflected in the strong resemblance of their 
observed spectra and hence they are normally investigated 
jointly. The observational data of this object were obtained 
from pe Marco et o/| ( |1997[ ) where 13 hues, given in Table [2j 
were extracted following a selection process. The plot 
indicates T ~ 16200 K which agrees very well with some 
previously-deduced values notably those of De Marco et al 
( De Marco et a/|1997||be Marco fc Barlow||2001 1. 



4 ELECTRON ENERGY DISTRIBUTIONS 



It has been suggested ( Nicholls et a^20"l2 l that the discrep- 
ancy between the results of ORLs and those of CELs is based 
on the assumption of a Maxwell-Boltzmann (M-B) for the 
electron distribution in the nebulae and that by assuming a 
different type of distribution, e.g. a K-distribution, the ORLs 
and CELs might yield very similar results for the abundance 
and electron temperature. One way for testing this proposal 
is to use DR lines to directly sample the distribution and 
compare to the M-B and other distributions. 

The cross-section for dielcctronic recombination can be 



expressed (e.g. Davies fc Seaton| ( 1969 1) 

_ IT / h 
"""■^ " 2 



(2) 



where P(e) is the dielectronic recombination probability, the 
probability of capture of a free electron of energy e and ve- 
locity V by an ion of statistical weight lu+ with the emission 
of a photon via a resonance state of statistical weight Ur; 
h and m being the reduced Planck's constant and mass of 
electron respectively. The recombination coefficient for di- 
electronic recombination via the resonance r is given by 



(3) 



where /(e) is the fraction of free electrons per unit energy. 
So 



2 CJ+ 



(-) fie,.) / Pie) de 



(4) 



where we have assumed that the resonance is narrow com- 
pared to changes in v and -P(e) so that Vr and e^ are the 
electron velocity and energy at the resonance position. Sev- 



eral workers have shown (e.g. Bell & Seaton ( 1985 1) that for 
an isolated narrow resonance 



Pie) de = 



2nh V T" 



(5) 



where F"^ and F" are probabilities of radiative decay and 
autoionisation in units of inverse time. We define a departure 
coefficient br by 



pr pa 



(6) 



which tends to unity for F" 3> F"^. We have assumed here 
that only dielectronic capture, autoionisation and radiative 



decay determine the population of the autoionising state. In 
principle there will also be radiative cascading from energet- 
ically higher autoionising states. In practice it is negligible 
for the transitions we consider here. The dielectronic recom- 
bination coefficient for a transition of wavelength A from the 
resonance state r is then given by 



1/2 



F'-(A) 6„ /(e,) 



(7) 



where ao is the Bohr radius, R is the Rydberg constant 
and where F''(A) is the radiative transition probability cor- 
responding to the line of wavelength A. The emissivity of a 
DR line is given by 

1 hp 

m = ^ N+ a^^iX) - (8) 

where A'''^ and A'"'" are the number density of electrons and 
ions respectively, h is the Planck's constant and c is the 
speed of light. Hence we may write for the fiux of the tran- 
sition with wavelength A 



7(A) = C 



,.(A) 



A 



where C is a proportionality factor and therefore 



fier) = DliX) 



1/2 A 



(9) 



(10) 



where D is another proportionality factor. The observed in- 
tensities are taken from Table [2] and the necessary atomic 
parameters are given in Table [l]Thus for each DR line aris- 
ing directly from an autoionising state (FF and FB tran- 
sitions), we can obtain the fraction of free electrons at the 
energy of that state. The resulting values of /(e) can then 
be compared to various theoretical electron energy distribu- 
tions. 

We derive values of /(e) for all the data sets that contain 
more than one FF or FB transition. In Figures [5|9| the results 
are presented. In these plots of e^2/(e) against e, the M-B 
distribution, 

2 



/mb (T, e) 



ikT) 



3/2 



e _ 
— e 

TV 



(11) 



appears as a straight line, where k is the Boltzmann con- 
stant, and T is the electron temperature shown at the opti- 
mum temperature obtained previously. 

We also show in these figures non-Maxwellian k- 
distributions for comparison, defined by ( |Bryans|200"5| 



(e) 



,3/2, 



F(k+ 1) 
I) 



F U- 



+ 

V «:e„ / 



(12) 



where k is a parameter characterising the distributions, 
while e^ is a characteristic energy. The K-distributions in 
these figures were calculated using the best fit temperatures 
from section 3. We show curves for k = 5 and k = 15, 
the former to illustrate clearly how the shape of a /t dis- 
tribution deviates from M-B in this representation, and the 
latter as representative of the values proposed by [Nicholls et\ 
\al (2012 \ to resolve the CEL/ORL abundance and tempera- 
ture problem. In the K-distribution there are more electrons 
than M-B at low energies and high energies and fewer at 
intermediate energies. The deviation of the K-distribution 
from M-B increases as k decreases. In typical nebular con- 
ditions the K-distribution is greater than M-B for energies 
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greater than 0.25-0.30 Ryd. This energy range is not ac- 
cessible with the lines and atomic data in use here. The low 
energy crossover occurs for energies typically below 0.05 Ryd 
in the Figures shown here. There are three groups of tran- 
sitions whose upper states lie in this energy region and can 
therefore potentially determine the shape of the distribution 
at the lowest energies and differentiate between the differ- 
ent distributions. They are 3d' ^P°-3p' lines 21-24 in 
Tablefl] 3d' 2D°/2-3p' ^Pg/j, line 38, and 3d' ^F°-3p' ^D", 
lines oT-53. The transition with the upper state with the 
lowest energy is line 38 at 6098.51 A which is only present 
in the spectrum of IC 418. Indeed all the spectra except 
IC 418 show only one multiplet in the low energy region, 
so that the distribution is poorly constrained. However, in- 
spection of Figures [5|7| suggests that, for the planetaries in 
our sample, the overall curvature of the K-distributions does 
not match the observations as well as the M-B distribution. 
In particular, in the case of IC 418, the object for which the 
observational data have the greatest number and spread of 
points, there appears to be significant depletion of electrons 
at the lowest and highest energies sampled relative to M-B 
and hence to any K-distribution. 

A more quantitative assessment of the K-distribution 
function that best matches the observational data can be 
obtained by calculating from the derived /(e) and the k- 
distribution as a function of k. However, the K-distribution 
is also a function of the characteristic energy e„ from which 
we can define a temperature from = e^jk which tends to 
the M-B temperature as k — >■ cxa. Hence the optimum match 
of a K-distribution may not occur at the temperature derived 
assuming a M-B distribution. Therefore we evaluate as 
a function of 1/k and r„ and the results are shown as con- 
tour plots in Figures 10-13. The errors on the observational 
data were derived as described in Section 2.3 above. Fig- 
ures [lO] and [Tl] show this measure for NGC 7009, the object 
with the largest ADF and CPD - 56°8032, an object where 
there is no reason to expect deviations from M-B. Note that 
in these plo ts 1/k = corresponds to a M-B distribution. 



4.1 Two-component models 



Figures 



12 



and 



13 



show the distributions for two more 
planetary nebulae with smaller ADFs than NGC 7009. For 
these latter two objects, the values of x'^ derived using the 
authors' flux error estimates were unreasonably large, so we 
have normalised x^ so that the minimum value as a function 
of K and is unity. We omit a contour plot for He 2-113 
since it is essentially the same as that for CPD - 56°8032. 
The two main features of these four figures are; 



(i) The minimum value of x^ occurs for 1/k = in all 
cases, so that the best value of — T. This is a quanti- 
tative reflection of the qualitative observation that the data 
points suggest negative curvature of f{e)/^ with respect 
to e whereas the M-B distribution has zero curvature and 
K-distributions have positive curvature. 

(ii) For values of close to the minimum and 1/k 7^ 0, 
differs little from its minimum value. Indeed, if we take 

a change of reduced x^ of unity to estimate the confldence 
interval on k, then 10 ^ k 5C 00 is possible in all cases. This 
is simply a reflection of the magnitude of the uncertainties, 
both in the observations with the current data sets extracted 
from the literature and in the calculated atomic parameters. 



Various authors (e.g. Liu et al (2000l) have suggested two- 



component models to explain the ORL/CEL abundance dis- 
crepancy, in which a cold metal rich component is embedded 
in a hotter medium. Suppose a nebula consists of two dis- 
tinct but individually homogeneous components (1 and 2) 
in which the electron energy distributions are assumed to 
be Maxwell-Boltzmann (/mb(2^i^))- The components have 
different electron number densities {Nf and A^l) and elec- 
tron temperatures (Ti and T2) and P is the fraction of the 
total volume {V) occupied by component 1. Then the power 
radiated in a line at wavelength A will be given by 



hr 

Pi(A) = F/3iVriV+a„„(A,ri)- 
P2(A) = V{1 - p)N^N+a^^{\,T2) 



he 
T 



(13) 



(14) 



where and A*'^ are the number densities of C^^ in com- 
ponents 1 and 2 respectively. Then adding Pi and P2 to get 
the total power, which is proportional to the observed flux, 
7, and using the expression for a^^ in terms of an electron 
energy distribution, we get 



7A 



-e^ (X /?Af A+/„3(e,Ti) + (l-/3)A|A+A„(e,T2)(15) 



where e is now the energy of the resonance which is the 
upper state of the transition of wavelength A. To plot the 
quantity on the RHS we need to relate the number densities 
in the two component model to those in a single component 
model for which 
JA 



(16) 



We require that if Ti = r2 = T the RHS are identical so 

(17) 



N^N+ = 13 NlN^ + (1 - ^) AIA2+ 



(18) 



where the a's are the abundances of C'^"'' relative to N"^ 
in the various components. The normalised two-component 
distribution is then 



/„3(e,r2) (19) 



The constants preceding the MB distribution functions are 
essentially relative fractional emission measures for the two 
components. 

We construct an illustrative two component model for 
NGC 7009 based on the model IH2 of |Liu et a/| p000| for 
NGC 6153, with Tl = 10000 K, T2 = 500 K and T = 5700 K 
from our work. We also assume Af = 5000 cm~^, A| = 
1000 cm-^ ai/a = 0.6, 02/0 = 70 and /3 = 0.7. The value 
of a2 = A(C^+)/A| is taken from the N{0^+)/N{R) ratio 



of Liu et al { 2000 1 in the metal-rich component of NGC 6153. 



With these parameters, the normalised two component dis- 
tribution is 



0.333 (e, Tl) + 0.667 /^^ (e, T2) 



(20) 



The resulting distribution is shown in Figure |5] Note that 
NGC 7009 has a significantly lower ADF (3-9, see Tab le [3l) 
than NGC 6153 (approximately 10, |Liu et al\ ( |2000[ )), so 
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Figure 5. Electron distribution plot for the NGC 7009 data 
of [Fang &: Liu| | |2011[ l, showing a Maxwell-Boltzmann and two 
K-distributions (re = 5.0 and 15.0) for T = 5700 K. A two- 
component Maxwell-Boltzmann distribution with Ti = 10000 K 
and T2 = 500 K is also plotted. The y-axis has an arbitrary scal- 
ing. The factor used to scale the observational data points to the 
theoretical curves was obtained by minimising the weighted least 
squares difference between the observational data points and their 
counterparts on a M-B distribution. 



departures from a single temperature model would be ex- 
pected be less extreme than modelled here. In addition, the 
emission measures and temperatures of the two components 
would be expected to differ between the two nebulae. Choos- 
ing a higher temperature (Ti) and larger emission measure 
for the low temperature component, for example, would im- 
prove the agreement with the low energy data point. Hence 
plots such as Figure |5] should place significant and useful 
constraints on the parameters of any two-component model. 



Figure 6. Electron distribution plot for NGC 7027 data of 
Zha ng et al\ | |2005[ l, showing a Maxwell-Boltzmann and two re- 
distributions (re = 5.0 and 15.0) for T = 12300 K. The ?/-axis 
scaling is as in Figure [5] 
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5 CONCLUSIONS 

We have determined an electron temperature from dielec- 
tronic recombination lines in a number of astronomical ob- 
jects, mainly planetary nebulae, by analysing C 11 transitions 
which originate from low-lying autoionising states and the 
subsequent decays. In the analysis we used a least squares 
minimisation method to find the electron temperature which 
best fits all the reliable observational data. For the plane- 
tary nebulae, our results generally fall below those derived 
from CELs such as the [O ill] forbidden lines and above 
those from ORLs like the O 11 permitted lines. There are 
exceptions however. In NGC 7027, for example, our temper- 
atures broadly confirm the forbidden line results while the 
O II ORLs yield a much lower temperature. In those objects, 
which are not planetary nebulae (DQ Her, CPD - 56°8032, 
and He 2-113), where a similar approach has been used be- 
fore, we find good agreement with the earlier results. We 
also find that the theoretical line emissivities that we pre- 
dict arc entirely consistent with those previously published 
by Davey et a/ (2000) for A4267. Given the very different but 
relatively simple mechanism of formation of the DR lines we 



10"! • • • • • • 1 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 

Electron Energy (Ryd) 

Figure 7. Electron distribution plot for IC 418 data of |Sharpee et\ 
[al| l |2003| , showing a Maxwell-Boltzmann and two K-distributions 
(re = 5.0 and 15.0) for T = 8200 K. The y-axis scaling is as in 
Figure |5] 

can conclude that the theoretical emissivities for A4267 and 
other similar C 11 lines are reliable. 

We have also proposed and demonstrated a method to 
test directly whether the free-electron energy distribution 
in planetary nebulae departs from Maxwell-Boltzmann. We 
showed that the fiuxes of DR lines originating directly from 
autoionising states can be used to sample the free-electron 
energy distribution and we applied this method to our sam- 
ple of objects. We showed that, for all the objects where 
suitable data are available, a Maxwell-Boltzmann distribu- 
tion gives the best fit to the observations but that the un- 
certainties in the observational data and atomic parame- 
ters are such that a K-distribution with values of k as sug- 
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Figure 8. Electron distribution plot for CPD - 56°8032 data of 
|De Marco et aZ| ( |1997[ l , showing a Maxwell-Boltzmann and two 
K-distributions (k = 5.0 and 15.0) for T = 17300 K. The y-axis 
scaling is as in Figure [5] 
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Figure 9. Electron distribution plot for He 2-113 data of Dc 
[Marco et al\ l |1997[ l, showing a Maxwell-Boltzmann and two re- 
distributions (k = 5.0 and 15.0) for T = 16200 K. The y-axis 
scaling is as in Figure [5] 



gested by Nicholls et al (20121 is not excluded. Similarly 
a two-component model as described by Liu et al ( 2000 1 



for NGC 6153 is also not excluded, although significant con- 
straints are imposed on any such model by the DR lines. We 
have highlighted several spectral lines which make it possi- 
ble to sample the low energy part of the electron energy 
distribution where departures from a simple one-component 
Maxwell-Boltzmann distribution are expected. Higher preci- 
sion is needed in the measurement of the intensities of these 
lines to fully realise the potential of this method. 
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Figure 10. derived from the difference between /(e) and a 
K-distribution as a function of 1/k and Tk for NGC 7009. 
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Figure 11. derived from the difference between /(e) and a 
K-distribution as a function of 1/k and for CPD - 56°8032. 
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Table 1: Observed transitions used in the current study where the columns stand for: an arbitrary transition index, the 
wavelength in A, the lower level designation and its statistical weight, the upper level designation and its statistical weight, 
transition status (FF, FB or BB), the radiative transition probability in s~^, the departure coefficient of the upper autoionising 
state in the case of FF and FB transitions and its energy in Ryd above the C^^ ls^2s^ ionisation threshold. The wavelengths 
are given as the vacuum's for A < 2000 A and as the air's for A > 2000 A. The prime in the level designation indicates an 
excited core, i.e. ls^2s2p(^P°), and the Is^ core is suppressed from all other configurations. 
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34 


5259.66 


3d' ^F° 


4 


4p' ''D'= 


2 


FB 


2.482E+7 


0.0167 


0.164896 


35 


5259.76 


3d' ^F° 


6 


4p' "D" 


4 


FB 


2.010E+7 


0.9429 


0.165014 


36 


5485.91 


3d' 


6 


4p' ^D" 


6 


FB 


3.323E+6 


0.9576 


0.165218 


37 


5648.07 


3s' *P° 


4 


3p' 


4 


BB 


1.945E+7 






38 


6098.51 


3p' ^P" 


4 


3d' ^D" 


6 


FB 


5.026E+7 


0.0261 


0.016144 


39 


6250.76 


3d' 


6 


4p P 


4 


FF 


2.733E+7 


0.8281 


0.161889 


40 


6779.94 


3s' ■'P" 


4 


3p' 


6 


BB 


2.497E+7 






41 


6780.59 


3s' ^P" 


2 


3p' ''D" 


4 


BB 


1.486E+7 






42 


6783.91 


3s' *P° 


6 


3p' 


8 


BB 


3.542E+7 






43 


6787.21 


3s' ''P° 


2 


3p' ^D" 


2 


BB 


2.946E+7 






44 


6791.47 


3s' *P° 


4 


3p' 


4 


BB 


1.875E+7 






45 


6798.10 


3s' -^P" 


4 


3p' ^D" 


2 


BB 


5.803E+6 






46 


6800.69 


3s' ^P" 


6 


3p' ^D^ 


6 


BB 


1.042E+7 






47 


6812.28 


3s' ^P" 


6 


3p' ^D'' 


4 


BB 


1.708E+6 






48 


7112.48 


3p' ■'D" 


2 


3d' *F° 


4 


BB 


3.321E+7 






49 


7113.04 


3p' ■'D" 


4 


3d' *F° 


6 


BB 


3.558E+7 






50 


7115.63 


3p' ■'D" 


6 


3d' ■'F" 


8 


BB 


4.062E+7 






51 


8794.08 


3p' ^D" 


6 


3d' ^F° 


8 


FB 


2.034E+7 


0.9982 


0.030241 


52 


8800.28 


3p' ^D" 


4 


3d' ^F° 


6 


FB 


1.897E+7 


0.9981 


0.029860 


53 


8826.55 


3p' ^D^ 


6 


3d' 2f° 


6 


FB 


1.405E+6 


0.9981 


0.029860 
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Table 2: The observed flux data for the investigated objects where the first column is the index as given in Table [l 



other columns give the investigated objects. For NGC 7027, the first column belongs to the data of Baluteau et al ( 1995 1 and 



while the 



the second to the data of Zhang et al (20051. For NGC 2867, the first value is related to knot 1 and the second to knot 2. 



The given flux is the value normalised to H/3 = 100 value except for the CPD - 56°8032 and He 2-113 where it is given as the 
absolute value in units of erg.s^^.cm^^ and in multiples of 10"^^. The symbol 'xx' indicates that the previous value of flux in 
that column is shared by the indicated lines in that multiplet. 



In. 

~ 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 



NGC 
7009 



NGC 
5315 



NGC 
7027 



NGC 
7027 



IC 
418 



NGC 
2867 



DQ 
Her 



CPD 
56°8032 



He 
2-113 



0.8795 

XX 
XX 



0.6559 

XX 
XX 



0.0021 
0.0021 



0.018 
0.059 

XX 



0.026 
0.575 

XX 
XX 

0.004 
0.026 
0.031 
0.019 
0.009 



0.0252 

XX 

0.0069 

XX 
XX 

0.5712 

XX 
XX 

0.0086 



0.0016 
0.0016 



0.0211 



28 


0.0088 








0.0044 


29 








0.013 


XX 


30 








0.013 




31 




0.0039 






0.0040 


32 








0.009 


0.0046 


33 








0.009 


0.0031 


34 










0.0032 


35 










XX 


36 


0.0004 










37 










0.0014 


38 










0.0011 


39 










0.0015 


40 






11.8 


0.034 


0.0109 


41 


0.0070 








0.0055 


42 






2.1 


0.004 


0.0022 


43 


0.0052 




3.7 


0.008 




44 




0.0068 


4.9 


0.012 


0.0066 


45 






0.7 






46 




0.0028 




0.009 


0.0050 


47 






0.5 


0.001 




48 






4.7 






49 






XX 




0.0052 


50 










0.0043 


51 


0.0320 




11.8 


0.015 




52 


0.0224 











270 

XX 
XX 



0.814(1.246) 

XX 
XX 



29 

XX 
XX 



4 


1940 


2 


2710 


6 


1670 


3 


3390 





1850 





1002 


1 


4370 





7141 





7183 





3571 





7183 





3571 


3 


5920 


1 


7850 





5906 





1973 


3 


3840 


1 


1300 


2 


8880 





9646 



0.045(0.079) 



53 



1.9260 
1.3470 
0.0943 



0.9342 
0.6533 
0.0457 
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Table 3: The range of electron temperature, from ORLs and a sample of CELs, in Kelvin, derived from different species and 
transitions, of the investigated astronomical objects as obtained from the literature where BD stands for Balmer Discontinuity, 
PD for Paschen Discontinuity, and RF for Radio Frequency. The value in the first row of each type of transition represents the 
minimum and the second is the maximum. Our results, as derived in the current paper, are given as averages in the additional 
three rows of C ll where the optimal value is in the first row, while the lower and upper limits of the confidence interval are in 



the second and third rows respectively. Our first value for NGC 7027 belongs to the data of Baluteau et al ( 1995 1 while the 



second belongs to the data of ^Zhang et al (2005 I. Similarly, the two values for NGC 2867 correspond to the first and second 
knots respectively. In the last row of the table, the range of the abundance discrepancy factor (ADF) values of some of the 
investigated objects as found in the literature is given. More details about the temperature data and their references can be 



found in Sochi (20121 



Object 


NGC 


NGC 


NGC 


IC 


NGC 


DQ 


CPD 


He 




7009 


5315 


7027 


418 


2867 


Her 


56°8032 


2-113 


H i(BD) 


7200 


8600 


8000 


>15000 


8950 


450 








8150 


8600 


12800 




8950 


1000 






H i(PD) 


5800 
5800 




8000 
8000 












He I 


5040 


10000 


8200 


9800 


10250 




20100 






8000 


10000 


10360 


9800 


10900 




20100 




C II 








9600 
9600 




700 
1450 


12800 
21700 


13600 
17000 


This work 


















Topt 


5650 


6950 


11100(12250) 


8200 


14300(16000) 


1600 


17300 


16200 




4830 


4620 


9930(8590) 


7250 


12120(13090) 


1460 


14020 


13630 


^max 


6240 


11990 


12660(23180) 


9630 


17600(26560) 


1740 


23680 


20740 


[Nil] 


10800 


9090 


12300 


8200 


8800 


2400 


11000 






11040 


10800 


12300 


9600 


11750 


2500 


11000 




II 


1600 
1600 


4350 
8100 


7100 
7100 












[O III] 


8350 


7800 


9260 


7000 


10520 










10380 


18500 


14850 


11200 


11850 








[N ii]+[0 III] 










11600 
11850 








RF 








6600 
23000 










ADF 


3.0-9.10" 


1.2-3.4" 


1.8^= 




1.49-1.77" 









° Obtained from|Liu et a/| ( |2004[ ); [Tsamis et d| ( |2004| ); [Liu et g^] ( (2006| ); [Barlow et a/| ( |2006| ); pu| ( |2006| ); |Wang fc Liu| ( |2007| ); 
Tsamis et a!{^20m. 



Obtained from Tlsamis et a«| (|2004|; [Garci'a-Rojas & Esteban [(|2007|. 
Obtained from ,Liu et al, (,2004^ 
Obtained from 



Garcia- Rojas et al (20091. 
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